O 
(N 

>>: 



A family of two- variable derivative polynomials for tangent and 

secant 

Shi-Mei Ma * 

Department of Information and Computing Science, Northeastern University at Qinhuangdao, 

Hebei 066004, China 



■ Abstract 

I In this paper we introduce a family of two- variable derivative polynomials for tangent and 

' secant. We study the generating functions for the coefficients of this family of polynomials. 

-4— > ' 

' In particular, we establish a connection between these generating functions and Eulcrian 

polynomials. 
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> ' 
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^ ■ 1 Introduction 

Tlj" , Throughout this paper, we define y = tan(2;) and z = sec{x). Denote by D the differential 

. operator d/dx. Thus D{y) = and D{z) = yz. An important tangent identity is given by 

X ■ 

\^ \ In 1995, Hoffman [9] considered two sequences of derivative polynomials defined respectively by 

D^{y) = Pn{y) a.nd D'^iz) = zQn{y) 

for n > 0. From the chain rule it follows that the polynomials Pn{u) satisfy Po(^) = and 
P„+i(n) = {l + u^)Pl^{u), and similarly Qo{u) = 1 and Qn+i{u) = {l + u^)Q'n{u) +uQn{u). The 
first few of the polynomials Pniu) are 

Pi{u) = l + u^, 
P2(u) = 2u + 2u^, 
P^[u) = 2 + 8n2 + 6n^. 

Various refinements of the polynomials Pniu) and Qn{u) have been pursued by several authors, 
see [a H El ini [H E] for instance. 
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Let Sn denote the symmetric group of all permutations of [n], where [n] = {1,2,... , re}. 
A permutation vr = 7r(l)7r(2) • • • 7r(re) G 5„ is alternating if 7r(l) > 7r(2) < •••7r(n). In other 
words, 7r{i) < 7r(i + 1) if i is even and 7r(i) > 7r(i + 1) if i is odd. Let En denote the number 
of alternating permutations in Sn- The number En is called an Euler number because Euler 
considered the numbers £^2n+i- There has been a huge literature on Euler numbers (see [T7] for 
details). In 1879, Andre [T] obtained that 



X 

re! 



y + z = ^En 

n=0 

Since the tangent is an odd function and the secant is an even function, we have 



2/ = IZ^2n+lT5-— TJ and Z = Y,E2n 



(2re + l)! ^ ""(2n)! 

n=0 ^ ' n=0 ^ ' 

For this reason the integers E2n+i are sometimes called the tangent numbers and the integers 
E2n are called the secant numbers. 

Let E{x) = y + z. Clearly, £'(0) = 1. It is easy to verify that 

2D{E{x)) = 1 + E'^{x). (1) 

Consider the derivative of ([T]) with respect to x, we have 

2'^D'^{E{x)) = 2E{x) + 2E^{x). (2) 

By the derivative of ^ with respect to x, we get 2^D'^{E{x)) = 2 + 8E'^{x) + GE'^^x). We now 
present a connection between E{x) and Pn{u) 

Proposition 1. For n> 0, we have 2"L'"(£'(x)) = Pn{E{x)). 

Proof. It suffices to consider the case re > 3. We proceed by induction on re. Assume that the 
statement is true for n = k. Then 

2k+iD''+\E{x)) = 2D{Pk{E{x))) 

= 2Pl{E{x))D{E{x)) 
= {l + E\x))P',{E{x)) 
= Pk+i{E{x)). 

Thus the statement is true for A; + 1, as desired. □ 
In |14j we write the derivative polynomials in terms of y and z as follows: 

L{n-1)/2J Ln/2J 

D^{y) = ^n,fcy"~^^"^-22'=+2 and L>"(z) = for n > 1. 

fc=0 k=0 

In particular, we observed that the coefficients Wn,k and ^ have simple combinatorial in- 
terpretations. The coefficient Wn k is the number of permutations in 5„ with k interior peaks, 
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where an interior peak of vr is an index 2 < z < n — 1 such that 7r(i — 1) < 7r(z) > Tr{i + 1). The 
coefficient ^ is the number of permutations in Sn with k left peaks, where a left peak of vr is 
either an interior peak or else the index 1 in the case 7r(l) > 7r(2) (see [6] for instance). 

In this paper we are concerned with a variation of the above definitions. The organization 
of this paper is as follows. In Section [21 we collect some notation, definitions and results that 
will be needed in the rest of the paper. In Section O we establish a connection between Eulerian 
numbers and the expansion of {Dy)^{y). In Section^ we establish a connection between Eulerian 
numbers of type B and the expansion of [Dy)'^{z). In Section [51 we study some polynomials 
related to {yDY{y) and {yDY{z). 

2 Preliminaries 

For a permutation vr € 5„, we define a descent to be a position i such that 7r(i) > 7r(i + 1). 
Denote by des (vr) the number of descents of vr. Let 

n 

ireSn k=l 

The polynomial An{x) is called an Eulerian polynomial, while A{n,k) is called an Eulerian 
number. The exponential generating function for An{x) is 

1 — X 

Aix,z) = l + Y,Mx)- = - j7^. (3) 

n>l 

The numbers A{n, k) satisfy the recurrence 

A{n + 1, fc) = kA{n, k) + {n - k + 2)A{n, k - 1) (4) 

with the initial conditions A{0,0) = 1 and A{0,k) = for /c > 1 (see [H A008292] for details). 
The first few of the Eulerian polynomials An{x) are 

^0(2;) = 1,^1 (x) = x,A2{x) = x + x^,A3{x) = x + ix^ + x^. 

It is well known that 



(5) 



An{x) — x"'~^^An ( ~ ) • 

\xj 

An explicit formula for A{n, k) is given as follows: 

A{n,k) = Y,{-ir(''^A{k-iT. 

Let Bn denote the set of signed permutations of ±[n] such that vr(— z) = — 7r(i) for all i, 
where ±[n] = {±1, ±2, . . . , ibn}. Let 



Bn{x) = Y,B{n,k)x'^ = 

fc=0 TrSBn 



^dcSB(7r) 
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where 

des B = \{i£ [n] : uj{i - 1) > + 1)}| 

with 7r(0) = 0. The polynomial Bn{x) is called an Eulerian polynomial of type B, while B{n, k) 
is called an Eulerian number of type B. Below are the polynomials Bn{x) for n < 3: 

Bo{x) = l,Bi{x) = l + x,B2{x) = l + 6x + x^,B3{x) = 1 + 23x + 23x^ + x^. 

The numbers B{n, k) satisfy the recurrence relation 

B{n + l,k) = {2k + l)B{n, k) + {2n - 2k + 3)B{n, k-1), (6) 

with the initial conditions B(0, 0) = 1 and B{0, k) = for k > 1. An explicit formula for B{n, k) 
is given as follows: 

B{n, k) = ^{-ly ^ ^) {2k -2i + 1)" 

for < A; < n (see [7] for details). 
For n > 0, we always assume that 

{DyT+\y) = {Dy){DyT{y) = D{y{DyT{y)), 

{Dyr+\z) = {Dy){Dyr{z) = D{y{Dyr{z)), 
{yDr+\y) = {yD){yDr{y) = yD{{yDT{y)), 
{yDT+\z) = {yD){yDT{z) = yD{{yDr{z)). 
Clearly, {DyY{y + z) = {DyY{y) + {DyY{z). For n > 1, we define 

2n 

{DyTiy + z)=Y, J(2n, k)y^^-^ z^+\ 

In Section [3] and Section [H we respectively obtain the following results: 

J(2n, 2k - 1) = 2"^(n, k) for l<k< n, 

and 

J{2n, 2k) = B{n, k) for 0<k< n. 

Let Jn{x) = Y.lloJ{'^n,k)x'' for n > 1. Then xJ„(x) = 2'^An{x'^) + xBn{x^). Therefore, 
from [14! Theorem 3], we have 

xJnix) = {l + xr+'An{x). (7) 

Using ([7|), we immediately obtain the following result. 
Proposition 2. For n > 1, we have 

n 

{DyTiv + z) = {y + zT^' Yl ^(^' 

k=l 
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3 On the expansion of {Dy)'^{y) 

For n > 1, we define 

n 

(Dyriy) = E{n, k)y^^-^^+^ (8) 

k=l 

Theorem 3. For 1 < k < n, we have E{n, k) = 2^A{n, k), where A{n, k) is an Eulerian number. 
Proof. Note that D{y'^) = 2yz'^. Then E{1, 1) = 2A{1, 1). Since 



,2n-2A:+l ^2k+2 



{Dyr+\y) = D{y{Dyr{y)) = 2YkE{n,k)y^''-"'+^z^''+2j2(.n-k+l)E{n,k)y' 

k=l k=l 

there follows 

E{n + l,k) = 2{kE{n, k) + {n - k + 2)E{n, k - 1)). (9) 
Comparison of Q with (HD gives the desired result. □ 

Let 

n 

Fn{y) = (Dyriy) = Y,F{n,k)y^^+\ (10) 

By Fn+i{y) = D{yFn{y)), we obtain 

= (1 + y^)Fn{y) + y{l + y^)F',{y) (11) 

with initial values -Fo(y) = U- We define Fn{y) = 2"a„(y) for n > 0. In the following we present 
an explicit formula for an{y). 

Theorem 4. For n>l, we have 

n 

an{y) = Y.A{n,k)y''^-^^^\l + y^)\ (12) 
fc=i 

Proof. Combining ^ and p^ . we get 

n 

Fn{y) = E{n, k)y^^-\l + y2)n-fc+l 

fc=l 

n / 2 \ ^ 

= 2"y-i(l+yy+i^yl(n,fc) . 

It follows from ([5]) that 

n 

= YMn,k)y'-~"^^Hi + yr- 

k=l 

This completes the proof. □ 



Let an{y) = Ylk=o"'(''^^^)y'^''~^^- Equating the coefficients of y'^^ ^'^^^ on both sides of ([T2 
we obtain 

a{n, n-k) = ^ A{n, i). 



It follows from (jlip that 

a{n +l,k) = {k + l)a(n, k) + ka{n, k - 1). (13) 

We define Wn{x) = Y2=o k)x^^^ . By ([l3]), we have 

Wn+i{x) = {x + x')Wl,{x), (14) 

with initial values VFo(x) = x. The first few terms of Wn{y) are given as follows: 

Wi{x) = x + x"^, 

W2{x) = x + 3x^ + 2x^, 

Wsix) = x + 7x'^ + 12x^ + 6x^, 

W^ix) = x + ISx^ + 50x2 _^ gQ^4 _^ 24x^ 

The triangular array {a(n, /c)}n>o,o<fe<n is called a Worpitzky triangle, and it has been exten- 
sively studied by many authors (see (TUl A028246]). 

In view of (jl4p . it is natural to consider the expansion of the operator ((x + x^)Z))". We 
define 

n 

{{x + x'^)D)'' = "^Gn,k{x){x + x'^)''D^ forn>l (15) 

k=l 

Applying the operator (x + x^)D on the left of (|15|) . we get 

Gn+l,kix) = + 2x)Gn,k{x) + ix + X^)D{Gn,k{x)) + Gn,k-l{x). (16) 

On the other hand, since 

D^{{x + x^)D) = {x + x^)D^^^ + k{l + 2x)D^ + k{k - l)D^-^. 
Applying the operator {x + x'^)D on the right of (|15p . we get 

Gn+i,k{x) = k{l + 2x)Gn,k{x) + k{k + l)(x + x^)Gn,k+iix) + Gn,k~iix). (17) 
Comparison of (fT6]) with pT|) gives D{Gn,k{x)) = k{k + l)G„^fc+i(3;). Thus 

G„,.(x) = ^!(^^^'"'(Gn,i(x)). 

Thus degG„^fc(2;) = n — k. For = 1, set Gn{x) = Gn,i{x). Then (fTHl) reduces to 

Gn+i{x) = (1 + 2x)G,(x) + (x + x2)D(Gn(x)) 



with initial values Gi{x) = 1. We define = G{n, k)x'^ ^ . It is easy to verify that 

G{n + l,k) = kG{n,k) + kG{n,k-l) (18) 

with 1) = 1. The Stirling numbers of the second kind, denoted by S{n, k), may be defined 
by the recurrence relation 

S{n + l,k) = kS{n,k) + S{n,k-l) (19) 

with the initial conditions 5(0, 0) = 1 and S{n, 0) = forn > 1 (see [Ml A008277 ] for details). 
Comparison of (jl8|) with (jl9p gives the following result. 



Proposition 5. For 1 < k < n, we have G{n, k) = k\S{n, k). 

4 On the expansion of {Dyy^{z) 

For n > 0, we define 

n 

(Z)y)"(z) = ^i7(n,A;)y2"-2fcz2fc+i_ 

k=0 

Theorem 6. For < k < n, we have H{n, k) = B{n, k), where B{n, k) is an Eulerian numbers 
of type B. 

Proof. Clearly, i^(0,0) = 1. Note that D{yz) = y^z + z^. Then H{1,0) = B{1,0) and H{1, 1) = 
B{1, 1). Note that 



{Dy){DyT{z) = ^(1 + 2k)H{n, A:)2/2"-2fc+2^2fc+i ^ Y^{2n - 2k + l)H{n, k)y 



2n-2fe^2fc+3 



A:=0 fc=0 

Thus we obtain 

H{n + 1, /c) = (1 + 2k)H{n, k) + {2n - 2k + 3)H{n, k-1). 

Hence H{n, k) satisfies the same recurrence and initial conditions as B[n, k), so they agree. □ 

Let (Dyriz) = zfn{y). Using {Dyr+\z) = D{yzfn{y)), we get 

fn+i{y) = (1 + 2y^)fniy) + y(l + y')/;(y) (20) 

with initial values /o(y) = 1- The first few terms of fn{y) are given as follows: 

/i(y) = 1 + 22/2, 

/2(y) = l + 8y2 + 8/, 

fsiy) = 1 + 26y2 + 72/ + 48^/^ 

/4(y) = 1 + 80y2 + 464y^ + 768y^ + 384y^ 
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Set fn{y) = El=of{n,k)y^''. By m, we obtain 

f{n + l,k) = {l + 2k)f{n,k) + 2kf{n,k-l) for < A; < n, (21) 
with initial conditions /(0, 0) = 1, /(O, k) = for k > 1. It should be noted [HI A145901] that 

(/(n,0),/(n,l),...,/(n,n)) 
is the /-vectors of the simplicial complexes dual to the permutohedra of type B. 

5 Polynomials related to {yD)^{y) and {yDY{z) 

For n> 1, we define 

n n 

{yDny) = Y,M{n,k)y^'^-h"+^ and {yDy\z) = Y,N{n,k)y''^z'^-^^+\ (22) 

k=l k=l 

Theorem 7. For 1 < k < n, we have 

M{n + l,k) = {2k - l)M{n, k) + (2n -2k + 4)M{n, k - 1), (23) 
N{n + l,k) = 2kN{n, k) + (2n - 2k + 3)N{n, k - 1). (24) 

Proof. Note that 

n n 

(yDJ^'+^y) = {yD){yD)''{y) = ^{2k-l)M{n,k)y'^^-^ z'^''-'^^+'^+^{2n-2k+2)M{n,k)y'^^+^ z^"" 

k=l k=l 

Thus we obtain (p3|) . Similarly, we get (p^ . □ 
From ()23p and (I24p . we immediately get a connection between M(n, A;) and iV(n, /c). 

Corollary 8. For 1 < k < n, we have M{n, k) = N{n, n — k + 1). 
Let 

n n 

Mn{x) = ^ M(n, k)x^ and iV„(x) = ^ iV(n, A;)x''. 

k=l k=l 

Then 

M„(x) = x'^+'Nn 0) . (25) 

Let 

n n 

Rn{y) = {yDYXy) = Y,R{n,k)y''^+^ and zTn{y) = (yD)"(z) = z J^T(n, fe)y2^ (26) 

fc=0 k=l 

Using (j26p . it is easy to verify that 

Rn+i{y) = y{i + y^)K{y) (27) 



and 



Tn+i{y) = y^Tniy) + y{l + y^^y). 



(28) 



The first few terms of Rniy) and Tn{y) are given as follows: 



Ri{y) 


1 3 

= y + y, 


J } ( ^ A 

R2Ky) 


= y + iy + 6y , 


RM 


= y + 13y3 + 27y^ + 15y^ 


Riiy) 


= y + 40y3 + 174y5 + 240y^ + 105y^ 


R^iy) 


= y + 121y2 + 990y^ + 2550y^ + 2625y^ + 945y 


Ti{y) 


= y^ 


T2{y) 


= 2y2 + 3/, 


T3iy) 


= 4y2 + 18/ + 15y^ 


n{y) 


= 8y2 + 84/ + 180/ + 105/ 


Uy) 


= 16y2 + 360/ + 1500/ + 2100/ + 945y^°. 



Equating the coefficient of ^y^fc+i both sides of (j27p . we get 

i?(n + l,k) = {2k + l)R{n, k) + {2k - l)R{n, k - 1). 

Equating the coefficient of y^'^ on both sides of (j28p . we get 

r(n + 1, k) = 2kT{n, k) + {2k - l)T(n, k - 1). 

Clearly, R{n, n) = T{n, n) = {2n — 1)!!, where {2n — 1)!! is the double factorial number. It should 
be noted that the triangular arrays {R{n, k)}n>ifl<k<n and {T{n, A;)}n>i,i<fc<n are both Gallon 
triangles (see [16l A187075 ] for instance), and it has been studied by Neuwirth [15]. We now 
present the following result. 

Theorem 9. For n > 1, we have 



Proof. Note that = + 1 . Combining (j22]) and ([26]) , we obtain 




and 




R^{y) = Y,M{n,k)y^^~\y^ + 1) 



k=l 



n-k+l _ 




and 




Using (j25]) . we get 




as desired. 



□ 
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By Theorem [91 we get i?„(l) = 7V„(2) and T„(l) = 2''iV„(i). It follows from that 

Nn+i{x) = (2n + l)xNn{x) + 2x(l - x)iV;(x) (29) 

with initial values Nq{x) = 1. The first few terms of Nn{x) can be computed directly as follows: 

Ni{x) = x, 

N2{x) = 2x + x^ 

N^{x) = 4x + IQx^ + x^, 

Ni{x) = 8x + mx^ + 36x^ + x^, 

A^5(x) = 16x + 296x2 + 516x^ + 116x^ + x^ 

In particular, N{n^ 1) = 2""^, N{n,n) = 1 and iV„(l) = (2n — 1)!! for n > 1. In the following 
discussion, we will study some properties of Nn{x). 

The numbers N[n, k) arise often in combinatorics and other branches of mathematics (see |12] 
for instance). A perfect matching of [2n] is a partition of [2n] into n blocks of size 2. Using (j24p 
and analyzing the placement of 2n — 1 and 2n, it is easy to verify that the number N{n, k) counts 
perfect matchings of [2n] with the restriction that only k matching pairs have odd smaller entries 
(see [m A185411]). It is weh known [HI A156919] that 



Ar„(x) = X;2--Q'), 



k\S{n,k)x''{l-x)''-'' for n > 1, (30) 

k=l ^ ' 

where S{n^ k) is the Stirling number of the second kind. By (jSOp . we get 

N{n, k) = ^(_i)'=-2"-2^ (^^ - .!5(n, i). 

Let 

4-n 

N{x,t) = TNn{x)- 

n>0 

Using (f29l) . the formal power series N{x,t) satisfies the following partial differential equation: 

^ .dN(x,t) , ^dN(x,t) , 
(1 - 2xt) - 2x(l - x) — = xN{x,t). 

By the method of characteristics [TB] , it is easy to derive an explicit form: 



iV(x,t)=e"* ' ^ ^ 



Hence 

Combining ([3]) and ([3T]) . we get the following result. 
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Theorem 10. For n > 0, we have 



(^V,(x)iV„_fc(x) = 2M„(x). 

L.—n V / 



k=0 

In the following of this section, we will give both central and local limit theorems for the 
coefficients of Nn{x). As an application of a result \13\ Theorem 2] on polynomials with only real 
zeros, the recurrence relation (I29p enables us to show that the polynomials {Nn.{x)}n>i form a 
Sturm sequence. 

Proposition 11. For n > 2, the polynomial Nn{x) has n distinct real zeros, separated by the 
zeros of Nn-i{x). 

Let {a(n, A;)}o<fc<n be a sequence of positive real numbers. It has no internal zeros if and only 
if there exist no indices i < j < k with a(n, i)a{n, k) ^ but a{n,j) = 0. Let An = X^fc=o "(^' 
We say the sequence {a{n,k)} satisfies a central limit theorem with mean and variance 
provided 



lim sup 

I— >+oo,a;Sl 



^ A „ , /t^ / 



An ^ 

The sequence satisfies a local limit theorem on i3 € M if 

cr„,a(n, [Xn + XOn) 1 



oo 



0. 



lim sup 



-e 2 



0. 



An 

Recall the following Bender's theorem. 

Theorem 12. 0/ Let {Pn}n>i be a sequence of polynomials with only real zeros. The sequence 
of the coefficients of Pn satisfies a central limit theorem with 



" Pnil) PuW \Pn{l). 

provided that lim = +oo. If the sequence of coefficients of each Pn{x) has no internal zeros, 
then the sequence of coefficients satisfies a local limit theorem. 

Combining Proposition [11] and Theorem 1121 we obtain the following result. 

Theorem 13. The sequence {N{n,k)}i<:k<n satisfies a central and a local limit theorem with 
fin = {2n + l)/4 and = (2n + l)/24, where n > 4. 

Proof. By differentiating ()29p . we obtain the recurrence Xn+i = {2n + 1)!! + (2n — l)xn for 
Xn = iV^(l), and this has the solution x„ = (2n + l)!!/4 for n > 2. By Theorem I12| we have 
^n = {2n + l)/4. Another differentiation leads to the recurrence 

(2n + l)!!, , , 

Vn+i = - 2) + (2n - 3)y„ 

for yn = N"{1). Set y„ = (2n — l)!!(an^ + 6n + c) and solve for a, b, c to get 

y„ = (2n - l)!!(12n2 - 8n - 7)/48 

for n > 4. Hence o"^ = (2n + l)/24. Thus lim a,^ = +oo as desired. □ 



n— )-oo 
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Let P{x) = ^i^QdiX^ be a polynomial. Let m be an index such that am = maxo<.j<n Oi- 
Darroch [5] showed that if P{x) E RZ(— oo,0], then 





< m < 




L^'n(l)J 









So the following result is immediate. 

Corollary 14. Let i = [{2n + l)J/4 or i = \{2n + 1)1/4. Then N{n,i) = maxi<fc<„ 7V(n, /c). 
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